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Abstract 

We investigate different choices for the quadratic momentum dependence required in nucleon- 
nucleon potentials to fit phase shifts in high partial-waves. In the Argonne v\g potential L 2 and 
(L • S) 2 operators are used to represent this dependence. The v\g potential is simple to use in many- 
body calculations since it has no quadratic momentum-dependent terms in S-waves. However, p 2 
rather than L 2 dependence occurs naturally in meson-exchange models of nuclear forces. We 
construct an alternate version of the Argonne potential, designated Argonne vis pq , in which the 
L 2 and (L • S) 2 operators are replaced by p 2 and Qij operators, respectively. The quadratic 
momentum-dependent terms are smaller in the Vis pq than in the v\s interaction. Results for the 
ground state binding energies of 3 H, 3 He, and 4 He, obtained with the variational Monte Carlo 
method, are presented for both the models with and without three-nucleon interactions. We find 
that the nuclear wave functions obtained with the vis pq are slightly larger than those with v±s at 
interparticle distances < 1 fm. The two models provide essentially the same binding in the light 
nuclei, although the vis pq gains less attraction when a fixed three-nucleon potential is added. 

PACS numbers: PACS numbers: 13.75.Cs, 21.30.-x, 21.45. +v 
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I. INTRODUCTION 



A number of choices have to be made in arriving at a form used to represent the nucleon- 
nucleon (AW) interaction by a model potential operator In general, models with different 
choices will give different results for systems with three and more nucleons. When we add 
many-body forces to the Hamiltonian, they too depend upon the choice made for Vij. In 
principle all realistic combinations of AA and three- nucleon (3 A) interaction, Vijk, should 
give identical results for the three-nucleon system by construction. We can similarly define 
combinations of 2-, 3-, 4- and many-nucleon interactions such that they give the same results 
for many-nucleon systems. However, such constructions are probably useful only when the 
contributions of n-body interactions decrease rapidly with n. In practice the contribution 
of realistic Vijk is much smaller than that of i>y, and the 4- and higher- nucleon terms in the 
Hamiltonian are often neglected. 

All fistic m odels of NN interaction HQS have a „„e-pion-e*ehan g e (OPE) character 
at long range, which is absolutely essential to fit high partial-wave data. The choices to 
be made for the OPE interaction are well known. As emphasized by Friar jj, 0], they are 
related by unitary transformations which generate the appropriate pion-exchange terms in 
Vijk- F° r example, a significant part of the difference between the Bonn potential Q and the 
Nijmegen ]J and Argonne j^J potentials is due to the choice made for the OPE interaction 

In contrast, the origins of the shorter-range parts of the AA interaction are not well 
known, and the choices made to model those are more pragmatic. When the model is 
intended for use in many-body problems, the accuracy with which it reproduces two-bod 
data and the ease with which it can be used in calculations are primary considerations 
Phenomenological three-body forces, determined from data are used with these models. 
They should reduce the dependence of the calculated many-body observables on the choices 
made to represent v^. 

Realistic models of the shorter-range parts of the AA interaction can be broadly cate- 
gorized as phenomenological, meson-exchange, or based on effective field theory. Here we 
focus on phenomenological models which can accurately reproduce the two-nucleon data. 
The meson-exchange models contain phenomenological effective mesons, and in models such 
as the CD Bonn [3J their masses and couplings depend upon the AA partial wave. Most 
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models based on effective field theory do not yet fit the scattering data as well. 
Phenomeno logical models use an operator construct, such as: 

V H = v p( r ij)°ij ■ (!) 

p=l,n 

They are based on the demonstration [8] that the most general form satisfying translational, 
rotational, parity, and time-reversal invariance can be written as: 

V NN = V c + V a <7i- (Tj + V t Sij + Vis (L ■ S)ij 

+v al L i:j + v ap (Ji • p a s ■ p , (2) 

with the tensor and quadratic spin-orbit operator definitions: 

= 3 Oi ■ f Oj ■ f — (jj ■ Oj (3) 
Lij = - fa -L,(Tj -L} . (4) 

Fortunately, on-energy-shell data require only two of the three possible tensor operators, 
, L^, and <Ji ■ p Oj ■ p . The last one has been omitted in all models Jjj] since the potential 

associated with it does not have a known theoretical motivation. 

The individual v x in Eq.Q can be functions of r, p 2 , or L 2 , and can have a general 

isospin dependence: 

v x = v° x (r, p 2 , L 2 ) + <(r, p 2 , L 2 )n ■ Tj (5) 
+ i;J(r,p 2 J L 2 )T ii +<-(r,p 2 > L 2 )(r zi + r zi ) , 

where the first two terms are charge- independent (CI), the next has the charge- dependent 
(CD) isotensor operator T^- = 3r zi r Z j — • Tj, and the last is a charge-symmetry-breaking 
(CSB) isovector component. The nuclear CD and CSB forces are much smaller than the CI 
terms. 

All modern AW potentials formulated in an operator structure contain the eight basic CI 
operators arising from the first line of Eq. (j2*j) : 

0£ =1 ' 8 = [1, a t ■ a„ S { j, L ■ S] ® [1, r t ■ Tj] . (6) 

These eight operators are sufficient to fit a CI average of S- and P-waves and to obtain good 
deuteron properties. Any interaction that is static or has linear terms in momentum can 



be represented as a sum of these eight operators. Typically six additional operator terms 
that are quadratic in either angular or linear momentum are required to fit higher partial 
waves. For example, the Urbana v u [l(3] and Argonne v u [ll] potentials use the additional 
operators: 

Of= 9 ' u = [L 2 , (a,, • er,)L 2 , (L • S) 2 ] ® [1, r, • r,-] , (7) 
where (L ■ S) 2 is simply related to Ly-: 

L tj = L • S - L 2 + 2(L • S) 2 . (8) 



The more modern Argonne v is (AV18) potential of Ref. |2| (hereafter referred to as WSS) 
supplements these fourteen CI operators with four small CD and CSB terms and a complete 
electromagnetic interaction. 

The S U pe,So ft -Co re (SSC) potent,. Q a,s„ used four L 2 operators, but with the an- 
gular momentum tensor operator Qij 

Qij = 3Lij — (cTj ■ £7j)L 2 , (9) 

or, in a later version [9], the Lij operator as the specific quadratic spin-orbit term. The Qij 
and (L • S) 2 operators vanish in S = states, while the does not. Still other variations 
of the quadratic spin-orbit operator, related via Eqs.flBJ) and Q, have been used in the past. 

The Urbana, Argonne, and SSC models are all local potentials in each NN partial wave 
denoted by quantum numbers L, S, J. In single channels having L = J, the v lsj is a central 
function of r, while in coupled channels having L = J ± 1, the vlsj is a sum of central, 
tensor and spin-orbit interactions. In contrast, the parametrized Paris potential was 
constructed with p 2 operators instead of L 2 , i.e., 

v ( r , p 2 ) = v c (r) + P V 2 (r) + v p2 (r)p 2 , (10) 

with 



P 2 = - 



r dr 2 r 2 



(11) 



where we have set % = 1 for brevity. Because p 2 does not commute with functions of r, a 
momentum-dependent potential like Paris can be more difficult to use in many-body cal- 
culations, especially configuration-space methods such as quantum Monte Carlo 14| and 
variational methods using chain summation techniques. In fact the Urbana model with L 2 
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and (L • S) 2 operators was proposed to avoid difficulties encountered in variational calcula- 
tions of nucleon matter with the Paris potential [15]. On the other hand, p 2 operators are 
more likely to arise in field-theoretic models for the NN interaction. 



In 1993 the Nijmegen group 



16| produced a multienergy partial-wave analysis of elastic 



NN data below Ti^ = 350 MeV (PWA93) that was able to reproduce over 4300 data points 
with a% 2 / datum ~ 1. This was followed by a series of comparably accurate potential models, 
including their own Nijm I, Nijm II, and Reid93 potentials Q, AV18 Q, and CD Bonn | . 
All these models use an operator construct in some fashion, with the Nijm II, Reid93, 
and AV18 choosing operators which give local potentials in each partial wave, while Nijm I 
includes p 2 terms and CD Bonn has additional momentum-dependent terms. However, with 
the exception of AV18, they are all adjusted partial-wave by partial-wave. The isoscalar t> 14 
part of the vis potentials contains only five independent interactions in S — 1 states and two 
in S = states, i.e., seven each for the two isospin values T = and 1. In the S = 1 states 
they correspond to the central, tensor, spin-orbit, L 2 or p 2 , and (L • S) 2 or Q^, while in the 
5 = states we have only the central and L 2 or p 2 potentials. In the T — 1, 5* = 1 states, 
for example, the five potentials can be extracted from the interactions in the 3 Po, 3 P\ and 
3 P 2 — 3 F 2 partial waves. The isoscalar interactions in all the 3 F 3 and higher T = 1, 5 = 1 
partial waves are related to those in the 3 Po, 3 -Pi and 3 P2 — 3 F 2 in a vi$ model. Similar 
relations exist in partial waves with S = 0. However, these relations are lost when the meson 
coupling constants or masses are varied to fit each partial wave, and have different values 
in each partial wave, as in the Nijmegen and Bonn models. These potentials therefore have 
additional nonlocalities which are best represented with partial-wave projection operators. 

Calculations of trinucleon and alpha binding energies with the five potentials mentioned 
above indicate that the three partial-wave local potentials give similar results, but the other 
two give more binding 0, Q|. CD-Bonn gives 0.4 and 2.0 MeV more binding than AV18 
in 3 H and 4 He, but a large part of it could be due to nonlocal representation of the OPE 
potential used in that model. As mentioned earlier, this extra binding goes into the three- 
nucleon interactions used with the other models. Nijm I uses the local OPE potential, but 
has shorter-range momentum-dependent terms. It gives only 0.1 and 0.7 MeV more binding. 
However, none of the models is able to reproduce the experimental binding energies without 
the addition of a three-nucleon potential. 

The purpose of this paper is to study the practical consequences of using p 2 instead of 
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L 2 operators in potential models. To do this, we construct a variant of the AV18 potential, 
using the same spatial functions and refitting the adjustable parameters to make the model 
as phase equivalent to AV18 as possible. As discussed below, our best reproduction of AV18 
phases also involves changing the quadratic spin-orbit operator from (L • S) 2 to Qy. As a 
result, we designate this new model as Argonne Vi 8pg (AV18pq) to denote the two changes 
in operator structure, L 2 — > p 2 and (L • S) 2 — > Qij. 

In Sec. II we describe the construction of the new potential and present a comparison of 
phase shifts and radial forms with respect to AV18. In Sec. Ill we present variational Monte 
Carlo calculations of the A = 3, 4 nuclei for both models, and with different three-nucleon 
potentials added. Our conclusions are given in Sec. IV. 

II. POTENTIAL CONSTRUCTION 

The practical need for the quadratic momentum- dependent operators is to simultaneously 
fit the phase shifts in S- and D-waves and in P- and F-waves. For example, a static potential 
with components of one- and two-pion exchange ranges and a shorter-range core cannot 
fit the 1 Sq and l D 2 phase shifts. The L 2 dependent term in the vis interaction operator 
generates a difference between the potentials in 1 Sq and 1 D 2 waves, and allows one to fit the 
phase shifts very accurately. 

Quadratic spin-orbit terms are necessary for fitting the variety of data in higher S = 1 
partial waves, such as 3 -Di,2,3 and 3 i 7 2,3,4- The AV18 is written with (L • S) 2 operators, 
but it can be algebraically recast with either the Ly or operators, with reshuffling 
of contributions among the (L • S) 2 , L • S and L 2 terms. As discussed below, we find it 
expedient to choose the operator, so we use AV18 in this representation, which we 
designate as AV18q for clarity. This will allow us to study the difference between choosing 
p 2 or L 2 keeping all other operators the same. This recast potential is the same as AV18, 
but its L 2 and L ■ S terms contain parts of the (L ■ S) 2 terms via Eqs. (JBJ) and (J5]). 

Inserting Eq.(JTTJ) in Ea.flTUl) we get 

v (r,p 2 ) = v c (r)-{--f^r,v p2 (r)} + 2v p2 (r)^ . (12) 

f (XT' f 

Therefore, in order to make comparison of terms in AV18q and AV18pq meaningful we 
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define: 

vfr(r) = \r*vf T (r) . (13) 

The L 2 term in AV18pq then has the same form as that in AV18q. Formally any v i$ pq model 
is obtained by adding the radial derivatives in p 2 , as given in Eq. ffTTj) . to a Vi 8q potential: 

= vis q - {~^r, irV 2 } . (14) 

However, the parameters of vi &q are not equal to those of Vi 8q which fits the phase shifts 
without the term with radial derivatives. 

The intermediate- and short-range parts of AV18pq are parameterized as in AV18 or 
equivalently AV18q. They are given by: 

v x s T (r) = vSr + r sT T*(r) (15) 



+ 



P x st + ^Q x st + ^YR x st W(r 



Here x can be c, t, Is, 12 or q for the central, tensor, spin-orbit, either L 2 or p 2 , and Qij 
interactions. The Vg T denotes the OPE interaction which contributes to terms with x = c 
and t. The average pion mass is /x, T M (r) is the tensor Yukawa function (with a cutoff), and 
the shape of the intermediate-range interaction is approximated with T 2 (r). A Woods-Saxon 
function W(r) is used to parametrize the short-range part. 

The radial wave equation for u(r)/r in single channels with the AV18q potential is 

where v c h{r) is the potential in that channel given by the Vi 8q operator. We solve this 
equation numerically using the Numerov method, which is particularly efficient for second- 
order differential equations without first derivative terms. With the addition of the derivative 
terms of Eq. (jl4|) . the radial wave equation for AV18pq generalizes to 

Q£ + 2 ^ T(r) K (r) + 2 ^(r)]V(r) (17) 

= (^^t^ + ^(r) - [vfAr)}" - E)u{r) . 

Although this expression has a first derivative of the wave function in it, it can be solved by 
the Numerov method with a simple substitution: 

w(r) = £(r)u(r) = [1 + AM T vf T (r)] 1/2 u{r) . (18) 



Then Eq. (jl7j) can be rewritten as 



D(r)w"(r) = r](r)w(r) , 



(19) 



with the definitions 



T](r) 




+ v ch (r) 



(jvfAr)}') 2 
D(r) 



— E . 



(21) 



(20) 



2M r r 2 



The above method is easily generalized for the coupled channels, because the new terms 
of the potential do not couple channels. In all other respects, the AV18pq potential is 
constructed with the same form as AV18q. The same set of fundamental constants is used. 
The long-range electromagnetic and one-pion-exchange terms are identical, including their 
short-range cutoffs. The strength parameters, J| T , Pf T , Q x STl and R X ST are refit to make 
AVI 8pq as phase equivalent as possible to AV18=AV18q. Note that the parameters -Pgr* 
and Qsr > are determined from the required behavior of v x (r — > 0), as discussed by WSS. 

The choice of Qy as the quadratic spin-orbit operator allows us to fit the S — and 
S = 1 phases for given T separately. The J < 4 phases are fit for eleven different energies, 
in the range 1 — 350 MeV, as well as scattering lengths and deuteron binding. In the T = 1 
phases, we want the quadratic momentum-dependent parts of the potential to be charge- 
independent, so we determine them first by fitting the phases produced by the CI part of 
the AV18 potential, including the CI average of one-pion exchange. Once these parts are 
fixed, we adjust the central part of the S = potential separately for pp and np phases, 
including full charge- dependence in OPE; this sets the CD parts of the interaction. The nn 
potential is then slightly altered from the pp potential to reproduce the nn scattering length, 
which is the only constraint on strong charge-symmetry breaking. The CSB in ST = 11 
channels is adjusted, as for AV18, such that the final potential in operator projection has 
no spin-dependent CSB, aside from electromagnetic terms. Fitting the T = channels is 
straightforward, including reproduction of the deuteron properties. 

The adjusted parameter values for AV18pq are given in Table HI in which lines labeled 12 
give the coefficients of the v l g T related to v P g T by Eq. ljlHj) . The rms deviations in degrees of 
the AV18pq phases from AV18 phases are given in Table HTl for pp scattering, and in Table ITTT1 
for np. These are binned over three different energy ranges, and in total. To help judge the 
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size of these errors, in the last column we quote the phase shift errors of PWA93 summed 
in a similar fashion. 

It is apparent from the table that the S = phases have been reproduced very well, with 
rms deviations much smaller than the experimental uncertainty. The 5 = 1 phases are more 
difficult to reproduce, although the rms deviations remain generally quite small. The places 
where they exceed the PWA93 errors are in the 3 P2, 3 F4, and higher T = 1 phases, and in 
the 3 Si, 3 D 3 , and 3 G4 T = phases. The deviations are concentrated at the higher energies. 

We also attempted fits using either (L • S) or L^j £ls alternate quadratic spin-orbit op- 
erators. Since (L • S) 2 vanishes in S = states, the potential in these channels is the same 
as in AV18pq. For the ST = 11 phases, we were able to get a slightly better reproduction 
of the AV18 model, including 3 P 2 and 3 F 4 deviations that fall below the PWA93 errors. 
However, the ST = 10 phases were notably worse, with 3 Di, 3 D2, and 3 G4 phases having 
total rms deviations ~ 0.3 degrees, which we judge to be unacceptable. To use Ly, we 
must fit S = 0, 1 channels simultaneously for given T, since this operator does not vanish 
in S — 0. We were able to get an even better reproduction of the AV18 T = phases with 
this operator, but we could not find acceptable T = 1 solutions. This does not mean that 
good p 2 potentials cannot be made with the or (L • S) 2 quadratic spin-orbit operators, 
but only that they may not work so well with the radial functions of the AV18 model. 

We now compare the AV18q and AV18pq interactions with each other and with the 
AV8' interaction which has no quadratic momentum-dependent terms. It has only the eight 
leading operators of Eq.fjHJ) plus Coulomb interaction between protons; it exactly reproduces 
the AV18 potential in 1 Sq, 3 Si — 3 D\ and all the P-waves, except for the small CD and CSB 
terms and non-Coulombic electroma gnet ic terms. The AV8' is obtained from the AV18 by a 
simple reprojection of the operators 19(, and is not phase equivalent to AV18 and AV18pq 
in D- and higher waves. Recall that AV18q is obtained by recasting AV18 with the identity: 

(L • S) 2 = \ (l + \a x ■ L 2 - ■ S + l -Q l0 , (22) 



and that the only difference between the operator structure of AV18q and AV18pq is L 



2 



P 2 . 



The central potentials in ST = 01 and 10 channels are compared in Fig. ^ In these 
spin-isospin channels the v c ST are those in the S and higher even waves, and they are the 
same in AV18q and AV8' models. The ST = 01 and 10 central potentials in AV18pq have 
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a slightly smaller repulsive core. This is to be expected. The L 2 and Qij operators are zero 
in S-waves. In the AV18q and AV8' models the change in the sign of S-wave phases at 
higher energies is produced by the repulsive core in the interaction. In contrast the p 2 part 
of AV18pq operates in the S'-waves also. It becomes more repulsive at larger energies, and 
helps turn the S-wave phase shifts negative together with the repulsive core. 

The Vg T in ST = 00 and 11 states are compared in Fig. they occur in P- and higher 
odd waves. In these states the AV8' has a larger repulsive core than either AV18 or AV18pq. 
The odd partial wave v c in AV8' contain the contribution of L 2 and Qij operators in P- 
waves. The repulsive cores in these states are obviously more sensitive to the treatment 
of quadratic momentum dependence as seen by comparing the AV18q and AV18pq curves. 
They are comparable to the centrifugal barrier, 2/ (Mr 2 ), in P- waves. For example, the 
centrifugal barrier is ~ 330 MeV at r = 0.5 fm. Fortunately nuclear energies are not very 
sensitive to the repulsive cores in P- and higher waves. 

The compared in Figs. El and |U Recall that the v p2 = r 2 i/ 2 /2, and that L 2 

also contains a r 2 factor. Therefore r 2 v l g T (r) are the relevant functions. In all cases the p 2 
dependent terms in AV18pq are smaller than the corresponding L 2 terms in AV18q. The 
Qij potentials are also smaller in AV18pq as shown in Fig. |S| 

The tensor and spin-orbit potentials are compared in Figs. |H] and The tensor potentials 
are weaker in AV8' but similar in AV18q and AV18pq. The large spin-orbit potential in 
ST =11 states is similar in all the models. However, the smaller spin-orbit potential in 
ST = 01 channel is relatively much weaker in AV18pq. 

The potentials in ST states can be cast into an operator format using <7j ■ <Jj and Tj ■ Tj 
operators as detailed in Sec. IV of WSS. The v c ST recast as the sum: 

v c (r) + v T (r)ri ■ tj + v a (r)(Ti ■ Oj + v aT (r)Ti ■ Tj a, ■ Oj (23) 

are shown in Fig. |H| At small r, where v c dominates, the smaller v T , v a and v ar are poorly 
determined by the data. The pair current operators depend upon these components of the 
interaction. For example, in Riska's method [2(| the isovector, pseudoscalar and vector 
currents are obtained from the v aT and Vt T potentials associated with the <7j • Gj 7j • Tj and 
Sij Ti ■ Tj operators. In particular the v a and v T are very different in the AV18pq and AV18q 
models. 

The deuteron properties of AV18q and AV18pq are compared in Table ITVl along with the 
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expectation values of different components of the potential. The binding energy and other 
observables are virtually identical. The total kinetic and potential energies obtained with 
AV18pq are smaller by ~ 1 %. The net contribution of quadratic momentum-dependent 
terms in AV18pq, (v p2 + v q ) = 0.086 MeV, is noticeably smaller than the (v l2 + v q ) = —0.331 
MeV for AV18q. 

The deuteron wave functions, u(r) and w(r), are illustrated in Fig. [HI the AV18pq func- 
tions peak at slightly smaller values of r. Impulse calculations of the deuteron structure 
function A(q 2 ) are virtually identical, but the B(q 2 ) structure function has its minimum 
pushed to larger momenta, and likewise the zero crossing of the tensor polarization T2o(q 2 ) 
also increases. However, consistent two-body current contributions constructed for AV18pq 
should be added before the electromagnetic observables are compared in detail with AV18 
predictions or experiment. 



III. LIGHT NUCLEI CALCULATIONS 



We have used variational Monte Carlo (VMC) techniques to calculate and compare the 
effect of AV8', AV18=AV18q, and AV18pq potentials on the binding of A = 3, 4 light nuclei. 



AV8' has been used in a benchmark study of four-nucleon binding |2l| and is the interaction 
used for propagation in Green's function Monte Carlo (GFMC) calculations. Finally, we 
study the effects of adding an explicit three-nucleon potential to the Hamiltonian, either 



Urbana IX (UIX) 



2J. (We use 



221 1 or the TM' variant |23J of the Tucson- Melbourne force 
TM' with the cutoff A = A.756m n as determined in Ref. [18] to be appropriate for use with 
AV18.) 

The VMC calculations are described in Refs. |23, bfj . We use the trial wave function 



l*v> 



i + E (u ijk + u™ 1 ) 

i<j<k 



l*P> 



(24) 



where the pair wave function, \l/p, is given by 



|#p> 



s]l(i + u i: 



i<j 



l*J> • 



(25) 



The Uij, Uijk, and U™ 1 are noncommuting two- and three-nucleon correlation operators, 
and S is a symmetrization operator. The antisymmetric Jastrow wave function, depends 
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on the state under investigation, but for S-shell nuclei takes the simple form: 



n fi 

i<j<k 



i]k\ r iki r jk) H Jc\' ijj 



Y[fc(T: 

i<j 



\$ A {JMTT 3 )) . 



(26) 



The fcijij) and f^ k are central (spin-isospin independent) two- and three-body correlation 
functions and $,4 is an antisymmetric spin-isospin state, 

1111. 



I $ 3(~~)) =A\ ]p[p]n) , 
|$ 4 (0000)) =A\]p{p]n{n) 



with A the antisymmetrization operator. 

The two-body correlation operator U{j is a sum of spin, isospin, and tensor terms: 



(27) 
(28) 



E 

P=2,6 



II fijk( T ik, r jk) 



uJr^O 



(29) 



where the Ofj are the first six operators of Eq.©. The f[j k is another central three-body 
correlation function. In previous calculations we also introduced a spin-orbit correlation in 
the sum of Eq.(J2U), but it is expensive to calculate and improves the binding energy only a 
little; we neglect such correlations here. 

The central fdr) and noncentral u p (r) pair correlation functions reflect the influence of 
the two-body potential at short distances, while satisfying asymptotic boundary conditions of 
cluster separability. Reasonable functions are generated by minimizing the two-body cluster 
energy of a somewhat modified NN interaction; this results in a set of coupled, Schrodinger- 
like, differential equations corresponding to linear combinations of the operators in v^, with 
a number of embedded variational parameters [2^. 

For the AV18pq potential, these differential equations must be generalized to take into 
account the influence of the v p2 terms. However, by studying the changes in the Schrodinger 
equation for the phase shifts, Eqs. (|16H2*T|) above, we find three simple substitutions that 
generalize the correlation generators, Eqs. (2.2-2.5) of Ref. j^j, to produce the appropriate 
correlation functions fsr'- 



2M r 
v c h(r) 
fsr(r) 



D(r) 



v c h(r - 



r\2 



D(r) 



fs T {r)/i{r) 
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A reprojection of the fsr(r) into operator terms then yields the / c (r) and u p (r). 

The ffj k , ffj k , and U^k are three-nucleon correlations also induced by details may 
be found in Ref. j^. The t/™ 7 are three-body correlations induced by the three-nucleon 
interaction, which have the form suggested by perturbation theory: 

U?jk 1 = Yl ' ' r v '• ('''./• fjk, hi) , (30) 

X 

with f = yr, y a scaling parameter, and e x a (small negative) strength parameter. 



Results of on]' VMC calculations arc given in Table V. We also show GFMG Q and 



Faddeev-Yakubovsky (FY) 18j results for comparison, where available. The VMC results 
are typically 1.5-3% above the more exact GFMC and FY results, which agree with each 
other to better than 1%. The primary goal here is to judge the relative binding of the 
different Hamiltonians, and we think the VMC calculations are adequate for this purpose. 

In all cases we see that the AV8', with no quadratic momentum dependence, gives the 
most binding. Both AV18 and AV18pq give less binding and the difference between their 
results is smaller than that from AV8'. From the point of view of fitting NN data, adding 
quadratic momentum-dependent terms to is the most important step, whether they be 
p 2 or L 2 . In A=3 and 4 nuclei these terms reduce the binding by ~ 2.5 and 4 % respectively. 
In the absence of three-nucleon interactions the reductions are nearly the same for AV18 
and AV18pq. (In examining the A = 3 results, it is important to remember that aside from 
Coulomb interaction, AV8' is a charge-independent potential, so averages of the 3 H and 3 He 
energies show the differences more clearly.) 

When a three-nucleon interaction is added to the Hamiltonian, as required to fit the 
experimental binding of light nuclei, the effect of adding quadratic momentum dependence 
grows to ~ 3 and 5 % for A=3 and 4 if AV18 is used, while a slightly larger effect of ~ 
4 and 6 % is obtained with AV18pq. In practice the strengths of the various terms in the 
three-nucleon interaction are not known a priori. The short-range terms of this interaction 
depend upon the representation of v^, and can be different for AV18 and AV18pq. They 
should be chosen separately such that both models give the same, experimental trinucleon 
energy. 

We note that this interplay between two- and three-nucleon interactions was already 
apparent in the Bochum FY calculations 3| for the Nijmegen potentials [|. These results 
are also shown at the bottom of Table IVl The nonlocal Nijm I potential gives slightly more 
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binding than the local Nijm II, but to obtain the experimental binding by adding a 3N 
force, a stronger version (larger A — ► weaker cutoff) is required to accompany the nonlocal 
NN potential, in this case chosen to reproduce the binding of 3 He. 

IV. CONCLUSIONS 

We have constructed an AV18pq two-nucleon interaction that is essentially phase equiv- 
alent to AV18. The difference between the two interaction operators is the choice made for 
quadratic momentum-dependent terms. The AV18pq uses the p 2 operator instead of the 
L 2 used in AV18. The results obtained with these two interaction models for the A = 2, 3, 
and 4 bound states are compared with those of the AV8' model which has no quadratic 
momentum dependence. 

The differences in the predicted deuteron wave functions are shown in Fig. E3 Those in 
the 4 He wave function can be judged from differences in the dominant central correlations 
in ST = 01, and central and tensor correlations in ST = 10 states shown in Fig. El In all 
these the magnitude of the wave function predicted by AV18pq is larger than that of AV18 
or AV8' at small interparticle distance r. This is due to the smaller repulsive cores in the 
S'-waves of AV18pq shown in Fig. ^ 

In principle the deuteron form factors predicted by AV18pq will be different from those 
of AV18. We plan to calculate the current operators for AV18pq following the methods used 
for AV18 Q, and compare the predicted form factors with available data. Such a comparison 
could indicate if one of the two choices is more realistic. 

The quadratic momentum-dependent terms give a small positive contribution to the 
ground-state energies of nuclei as can be seen from Table |3 This contribution can depend 
upon the choice of interaction operators. In the absence of a 3N interaction, the quadratic 
momentum-dependent terms in AV18 and AV18pq give almost identical corrections, within 
0.5%. However, a larger difference between AV18 and AV18pq develops when a fixed realistic 
3N interaction is added to the nuclear Hamiltonian to obtain binding energies closer to the 
experimental values. In this case, the non-local AV18pq becomes less attractive than the 
local AV18, by ~ 1%, just as the non-local Nijm I became relatively less attractive compared 
to the local Nijm II potential. One can presumably choose the parameters of the Vijk to 
reproduce binding energies of nuclei with either AV18 or AV18pq. 
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All one-boson exchange models of nuclear forces give interactions with a p 2 dependence. 
Electromagnetic and color forces also have a p 2 dependence. On the other hand, nucleons 
are composite objects with internal excitations, and their interactions can have either L 2 or 
p 2 terms or both. The main motivation behind the choice of L 2 was to simplify many-body 
calculations [15J. Green's function Monte Carlo methods can be used to study nuclei having 
A < 12 with AV8' type interactions that contain the dominant features of nuclear forces, 
but not the quadratic momentum dependence. 

In GFMC studies with the AV18 interaction, the quadratic momentum-dependent terms 
are treated as a first-order perturbation The error 5E in this calculation is given by: 

5E = (ttisliZislttw) - (* 8 ,|F 18 |* 8 ,) , (31) 

where | \I/ 18 ) and |^8') are the eigenstates obtained with the Hamiltonian if 18 having the 
AV18 interaction and with H 8 * having the AV8' respectively. The Hamiltonians can include 
any static 3N interaction. The values obtained for 5E in 4 He, using optimum variational 
wave functions to approximate the \^is) and 1^8') are respectively 0.05, 0.16 and 0.03 MeV 
for AV18 alone and with UIX and TM' models of 3N interaction. They are less than 1% of 
the total energy. 

A similar study using AV18pq instead of AV18 gives much larger estimates of 1.6, 1.1 and 
1.8 MeV for the SE in the above three cases. These larger values are due to differences in 
the wave functions of AV8' and AV18pq models at small r (see Figs. I§1 and irU|) . A first-order 
treatment of the difference between AV18pq and AV8' gives much too large an effect. It will 
be necessary to develop new techniques to perform accurate GFMC many-body calculations 
with the AV18pq interaction. 

For those who wish to use the AV18pq potential model in their own applications, a 
FORTRAN subroutine is available online \27\. 
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TABLE I: Short-range potential parameters in MeV. The asterisk denotes that the value was 
computed by Eq. (24) of WSS and not fit. The three shape parameters are: c = 2.1 fm -2 , vq = 0.5 
fm, and a = 0.2 fm. 



Channel 


Type 


i 


P 


Q 


R 


S = 0,T = 


l(pp) 


c 


-10.518030 


2836.0715 


1582.7028* 


651.1945 


S = 0, T = 


lffio) 


c 


-10.812190 


2816.4190 


1578.2721* 


1002.5300 


S = 0,T = 


1 inn) 


c 


-10.518030 


2832.4903 


1580.7610* 


651.1945 


S = 0,T 


= 1 


/2 


0.134747 


-9.4691 


-5.1342* 





S = 0,T 


= 


c 


-4.739629 


1121.2225 


466.8222* 


2764.3395 






/2 


-0.227084 


166.5629 


90.3117* 





5 = 1,T = 


l(pp) 


c 


-9.882847 


2589.7742 


1389.2081* 


2952.3910 


5 = 1,T = 


l(nv) 


c 


-9.882847 


2587.9836 


1386.1622* 


2952.3910 


5 = 1,T = 


l(?T.n.) 


c 


-9.882847 


2586.1930 


1387.2664* 


2952.3910 


5 = 1,T 


= 1 


/2 


-0.008159 


132.7694 


71.9886* 


-169.8510 








1.420069 





-453.5357 


-837.3820 






Is 


-1.749197 


-493.8470 


-267.7677* 


1533.0637 






Q 


0.135181 


-17.7975 


-9.6499* 


-46.2542 


5= 1,T 


= 


c 


-8.351808 


2325.5929 


1307.9923* 


957.8091 






12 


-0.023577 


1.8164 


0.9849* 


127.2921 






t 


1.327862 





-1170.8528 


580.5596 






Is 


0.060223 


58.3208 


31.6220* 


-126.0235 






q 


0.000589 


-25.1123 


-13.6161* 


-4.6897 
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TABLE II: RMS deviations 
Channel 



of pp phases 
1-25 50-150 
MeV MeV 



in degrees for AVT8pq compared to AV18. 
200-350 Total PWA93 
MeV 



1 a 


0.021 


0.031 


1 7-1 


0.001 


0.011 


1 /~i 





0.001 


3 r> 

*Po 


0.024 


0.024 


3 n 


0.002 


0.035 


3 p 


0.008 


0.090 


£2 


0.002 


0.011 


3 ^2 


0.001 


0.013 


3 ^3 


0.002 


0.042 




0.002 


0.054 


e 4 





0.005 


3 #4 





0.018 



0.245 0.149 0.181 

0.061 0.037 0.055 

0.020 0.012 0.029 

0.092 0.059 0.224 

0.074 0.048 0.109 

0.160 0.107 0.062 

0.044 0.027 0.048 

0.055 0.034 0.062 

0.034 0.030 0.054 

0.077 0.054 0.043 

0.009 0.006 0.001 

0.047 0.030 



TABLE III: RMS deviations of np phases in degrees for AV18pq compared to AV18. 
Channel 1-25 50-150 200-350 Total PWA93 





MeV 


MeV 


MeV 






'So 


0.018 


0.087 


0.078 


0.067 


0.344 


'Pi 


0.005 


0.029 


0.064 


0.041 


0.160 


l F 3 





0.003 


0.022 


0.013 


0.047 


3 s 1 


0.095 


0.113 


0.387 


0.248 


0.154 




0.011 


0.061 


0.048 


0.043 


0.115 


3 D, 


0.002 


0.086 


0.200 


0.129 


0.127 


3 D 2 





0.037 


0.139 


0.086 


0.090 


3 D 3 


0.001 


0.017 


0.230 


0.139 


0.069 


£3 


0.001 


0.005 


0.036 


0.022 


0.035 


3 G 3 





0.001 


0.004 


0.002 


0.002 


3 G 4 





0.027 


0.253 


0.153 
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TABLE IV: Comparison of static deuteron properties and expectation values for different terms in 
the potentials. 





Expt 


AV18q 


AV18pq 




E d 


2.224575(9) 


2.224575 


2.224573 


MeV 


(T) 




19.810 


19.525 


MeV 


(V) 




-22.035 


-21.750 


MeV 


(v c ) 




-4.446 


-4.079 


MeV 


(v*) 




-16.586 


-17.462 


MeV 


(v ts ) 




-0.690 


-0.313 


MeV 


1 ]1\ 

(v ) 




1.087 




MeV 


i r>9\ 

(v pZ ) 






1.440 


MeV 


(vi) 




— 1.418 


-1.354 


MeV 






0.018 


0.018 


MeV 


A s 


0.8781(44) 


0.8850 


0.8863 


f m l/2 


V 


0.0256(4) 


0.0250 


0.0250 




rd 


1.953(3) 


1.967 


1.969 


fm 


Vd 


0.857406(1) 


0.847 


0.847 


Mo 


Qd 


0.2859(3) 


0.270 


0.270 


fm 2 


Pd 




5.76 


5.78 


% 
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TABLE V: Binding energies of light nuclei (in MeV) for different Hamiltonians as computed by 
the variational Monte Carlo (VMC) method, with comparisons to relevant Green's function Monte 
Carlo (GFMC) results of Ref. Q and Faddeev-Yakubovsky (FY) results of Ref. Q- We also show 
FY results for the Nijmegen potential at the bottom. Different versions of the Tucson-Melbourne 
(TM) force are characterized by the cutoff parameter A shown in parentheses. 



Hamiitoman 


Method 


3tt 


3tt 

J He 


4tt 

He 


AV18 


VMC 


-7.50(1) 


-6.77(1) 


-23.70(2) 




GFMC 


-7.61(1) 


-6.87(1) 


-24.07(4) 




FY 


-7.623 


-6.924 


-24.28 


AV18pq 


VMC 


-7.50(1) 


-6.77(1) 


-23.79(2) 


AV8' 


VMC 


-7.65(1) 


-7.01(1) 


-24.69(2) 




GFMC 


-7.76(1) 


-7.12(1) 


-25.14(2) 


AV18/UIX 


VMC 


-8.29(1) 


-7.53(1) 


-27.58(2) 




GFMC 


-8.46(1) 


-7.71(1) 


-28.33(2) 




FY 


-8.478 


-7.760 


-28.50 


AV18pq/UIX 


VMC 


-8.22(1) 


-7.47(1) 


-27.21(2) 


AV8'/UIX 


VMC 


-8.51(1) 


-7.86(1) 


-29.05(2) 




GFMC 


-8.68(1) 


-8.03(1) 


-29.82(2) 


AV18/TM'(4.756) 


VMC 


-8.26(1) 


-7.50(1) 


-27.51(2) 




FY 


-8.444 


-7.728 


-28.36 


AV18pq/TM'(4.756) 


VMC 


-8.22(1) 


-7.46(1) 


-27.42(2) 


AV8'/TM'(4.756) 


VMC 


-8.44(1) 


-7.79(1) 


-28.83(2) 


Nijm I 


FY 


-7.736 


-7.085 


-24.98 


Nijm II 


FY 


-7.654 


-7.012 


-24.56 


Nijm I/TM(5.035) 


FY 


-8.392 


-7.720 


-28.60 


Nijm II/TM(4.975) 


FY 


-8.386 


-7.720 


-28.54 
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FIG. 1: (Color online) Central potentials in even-parity waves. 
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FIG. 2: (Color online) Central potentials in odd-parity waves. 
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FIG. 3: (Color online) Quadratic momentum-dependent potentials in even-parity waves. 
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FIG. 5: (Color online) Quadratic spin-orbit potentials in S = 1 channels. 
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FIG. 6: (Color online) Tensor potentials in S = 1 channels. 
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FIG. 9: (Color online) Deuteron S- and D-wave functions. 
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